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We derive thermodynamic Bethe ansatz equations for a model describing an U —* oo Anderson impurity 
embedded in a BCS superconductor. The equations are solved analytically in the zero-temperature limit, 
T = 0. It is shown that the impurities depress superconductivity in the Kondo limit, however at T = the 
system remains in the superconducting state for any impurity concentration. In the mixed-valence regime, 
an impurity contribution to the density of states of the system near the Fermi level overcompensates a 
Cooper pairs weakening, and superconductivity is enhanced. 
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Since the pioneering work of Abrikosov and Gor'kov 
(AG) [|J , the problem of superconducting magnetic alloys 
has been the subject of many early Q and more recent ||] 
studies. Almost all of the theoretical methods developed 
to attack the Kondo problem in normal metals, from per- 
turbative approaches to Wilson's numerical renormaliza- 
tion group, have been generalized to the case of supercon- 
ductors. Perhaps the only exception is the Bethe ansatz 
(BA) technique, which solves the Kondo problem in nor- 
mal metals exactly, but cannot be straightforwardly gen- 
eralized to the case of dilute superconducting alloys. The 
basic theoretical models describing magnetic impurities 
in normal metals, such as the s-d (Kondo) and Anderson 
models, are integrable under two additional conditions: 
(i) an electron-impurity coupling is assumed to be energy 
independent, and (ii) a band electron dispersion, Ej~, can 
be linearized around the Fermi level, Ek — vp(k — kp), 
where kp and vp are the Fermi momentum and veloc- 
ity, respectively Q| . Since a carrier dispersion in the su- 
perconducting state cannot be linearized near the Fermi 
level, these conditions eliminate superconductivity from 
the BA analysis of the behavior of magnetic alloys. 

However, it has recently been discerned || that the 
basic "impurity" models of quantum optics, describing a 
system of Bose particles with a nonlinear dispersion cou- 
pled to two-level atoms, exhibit hidden integrability and 
are thus exactly diagonalized by BA. One of the most ex- 
citing potential applications of the approach developed 
may be an extension of the BA method to the Kondo 
problem in superconductors and other Fermi systems (e. 
g., gapless Fermi systems |Q) with an essentially nonlin- 
ear dispersion of charge carriers. 

In this Letter, we employ hidden integrability of a 
model describing an Anderson impurity with an infinitely 
large Coulomb repulsion on an impurity orbital embed- 
ded in a BCS superconductor to study the thermody- 
namic properties of the system. In the standard manner 
, we find a set of basic thermodynamic equations for 
the energies of elementary excitations of the system. In 
terms of these equations, we derive an exact equation for 
the order parameter of the superconducting phase tran- 



sition |S|], A, minimizing the thermodynamic potential of 
the system, f2, with respect to A, Stl/SA = 0. While at 
finite temperatures the basic equations require a numer- 
ical analysis, in the zero-temperature limit, T = 0, they 
are solved analytically, giving an exact expression for the 
impurity contribution to the parameter A. 

At T = 0, the order parameter is given by the expres- 
sion A = Aoexp (— /Ximp), where Ao is the order param- 
eter in the absence of impurities and the parameter /ii mp 
[see Eq. (15)] describes the impurity contribution. The 
magnitude and sign of /ii mp are determined completely 
by the position of the impurity energy level, e^, with 
respect to the Fermi energy of the host metal, ep. In 
the Kondo limit, where lies much below ep, /x; mp is 
positive, and hence magnetic impurities depress super- 
conductivity. However, the system remains in the super- 
conducting state at any concentration of impurities. In 
the mixed-valence regime, where lies near the Fermi 
level, the parameter /Xi mp becomes negative. An impurity 
contribution to the density of states of the system near 
the Fermi level dominates over a Cooper pairs weakening, 
and the Anderson impurities enhance superconductivity. 

Throughout the paper, the BA technique and many 
results of the exact solution of the Anderson model in 
normal metals are often used with no special references. 
All needed details can be found in excellent comprehen- 
sive reviews jj] and references therein. 

We start with a Hamiltonian that includes the Hamil- 
tonians of the BCS and Anderson models, 

H = Yl E k a la a k<y ~ ( Aa U a kl + A*a fci a fcT ) 

k.a k 

+ A Vg + Y Vk (4<A + 4«fc CT ) 

k,cr 

+e d J24d, + Ud\d^d\d l , (1) 

a 

where we have used the standard spherical harmonic rep- 
resentation for band electron operators. The Fermi oper- 
ator at creates a conduction electron with the momen- 
turn modulus k, spin a =|, | and the energy Ek — e^ — ep, 
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where ek is the kinetic energy. Only the s-wave is as- 
sumed to be coupled to the impurity, therefore all other 
partial waves have been dropped. An electron localized 
in the impurity orbital is described by the Fermi opera- 
tors d a . The fourth term of Eq. (1) represents the hy- 
bridization of the band and impurity level electrons with 
the matrix element Vk, while the Coulomb repulsion on 
the impurity orbital is described by the last term. The 
parameter A = gX)fc( a fcJ. a fct) i s assumed to result from 
the Cooper pairing phenomenon with positive coupling 
constant g. 

Diagonalization of the BCS part of Eq. (1) by the 
Bogoliubov-Valatin unitary transform |1 gives 



H = E BCS + ^2 ^Aa c k<y + V Y1 (4 c fc^ + 4 CT d CT 



ka ko 

-e d ^2 d\d a + Ud\d^d[di, 



(2) 



where ui% — — \/k 2 + A 2 , for k < 0, — Vk 2 + A 2 for 
k > 0, and .Ebcs = J2k(^~ LU k) + A 2 /g- For simplicity, we 
have linearized the band spectrum of the host metal in 
the normal state around the Fermi level and set vp = 1. 
The electron momentum and energy are taken relative 
to the Fermi values. We have also set k = kp in both 
the hybridization matrix element, v — v(kp), and in the 
coefficients of the unitary transform. 

Moreover, to apply the BA method to the Hamiltonian 
(2), we have omitted the terms d^ct and d a Cka, which 
do not conserve the number of excitations in the system; 
these terms are assumed to lead only to insignificant cor- 
rections. The bare vacuum of the model is then defined 
by Cka\0) = d a \0) = 0. To obtain the ground state of the 
system in the absence of the impurity, one thus needs 
to fill all states with k < 0. In the normal state, A = 0, 
Eq. (2) reduces to the integrable version of the Anderson 
model diagonalized by Wiegmann |lC(| . 

In what follows, we confine ourselves to very large val- 
ues of U, e d + U > 2?(+) , where 2?(+) is the upper edge 
of the band, so that double occupancy of the impurity 
level is excluded. The eigenvalues of the model Hamilto- 
nian (2) are then found from the following Bethe ansatz 
equations (BAE) B: 



exp (ikjL) 



N 

n 

i=i 



hj -e d /2T-i/2 



hj - e d /2T + i/2 



M u 

a—1 J 



X a - i/2 



X a + i/2 



i/2 



■i/2 



M 

n 

/3=1 



X n — Xr 



X a — Xp + : 



(3a) 



(3b) 



where E = J^j^i^j) i s the eigenenergy, N is the to- 
tal number of particles in the interval L, and M is 
the number of particles with spin down. The function 

hj = h(kj) is defined by h(k) = + (l - ^j^y^ , 

where T — npv 2 — v 2 /2 and p = l/2n is the density of 
states of the host metal in the normal state. 



In the thermodynamic limit, spin "rapidities" X a are 
grouped into bound spin complexes of size n, 



A^" j) = X™ +i(n+l- 2j)/2, j = l,. 



(4) 



The simplest bound states of charge excitations are as- 
sociated with real spin rapidities X a , and their charge 
rapidities fci ± ' ) are found from the equation 



h(k^)=X a ±i/2. 



(5) 



In the normal state, Eq. (5) has a single solution, k^ = 
2r(A Q ± i/2). Since A T, this solution acquires in the 
superconducting state a small correction of the order of 

(A/r) 2 , 



fc( ± )(A)c2r(A± l /2)-i(^Ay 



(x±i/2y 



(6a) 



The energy of such "normal" charge complexes, o;„(A) = 
w(/c' + ' ) ) +uj(k^), is found to be 



^n(A) 

4r 



A + 



A - e d /2T 
A 2 + 1/4 ' 



(6b) 



where the second term describes the "gap" correction. 
The sign of this correction is different for different A, 
therefore the appearance of the energy gap can either 
increase or decrease the total energy of normal charge 
complexes. In what follows, this fact will play a crucial 
role in the interplay between the magnetic and supercon- 
ducting properties of the system. 

Eq. (5) also admits "gap" charge complexes with ra- 
pidities p a ^ — zkiAcos z a ^ and the energy w g (A Q ) = 
2 ARe sin z a . Here 



(Aa) 



*(-) 



-(A a ) = arctan ~^ 2T 
y a> A„ /or J 



, a -e d /2Y + i/2 



(7) 



and the terms of order A/r are omitted. In what fol- 
lows, we consider the case of negative e d . The gap com- 
plexes can then be shown to exist only for A > e d /2T, 
and their energies are positive, uj g > 0. BAE admit also 
"long" charge complexes associated with spin complexes 
(4). We however do not consider such excitations, be- 
cause, as in the normal state, they do not contribute 
to the low-temperature thermodynamics of the system. 
Finally, the subgap spectrum of the model contains a 
discrete mode (DM) ||, which is naturally treated as 
a particle-impurity bound state. In the BA approach, a 
DM is first found as a solution of a single particle problem 
rather than as a multiparticle discrete mode predicted by 
Shiba . It can be shown that for e d < the renormal- 
ized energy of a DM in the multiparticle spectrum of the 
model (2) is much bigger than A, and hence a DM does 
not contribute to thermodynamics of the system in the 
temperature range of physical interest, T < T c , where T c 
is the superconducting critical temperature. 
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In the standard manner, we then find in the thermody- 
namic limit a set of equations for the renormalizcd ener- 
gies of the elementary excitations e(k), £(A) and K n (X), 
corresponding to unpaired charge excitations with real 
kj, normal charge complexes and spin complexes, respec- 
tively: 



+ 



e{k) = Lu{k)+ / d\ ai (h(k)-\)F[-t(\)] 

' — CO 

d\ai(h(k) — A) F[— 77(A)] 
f d\a n (\-h(k))F[- Kn (X)} (8a) 

J — OO 



d /2r 



E 

n=l • 



£(A)=w„(A)+ / dkh'(k) ai (X-h(k))F[-e(k)] 

J — OO 

/CO 
dAa 2 (A- A')F[-<£(A')] 
-CO 
/>CO 

+ / dAa 2 (A- A')F[-r ? (A')]) (8b) 

°° />CO 

F[« n (A)] = J2 dXA nm (X-X')F[- Km (X')] 

/oo 
da'(fc) fln (A-/i(fc))F[-e(fc)]. (8c) 
-CO 

The renormalized energy of gap complexes is given by 

77(A) =w g (A)+£(A)-u;„(A), X>e d /2T. (8d) 

Here, F[f(x)} = Tin [1 + exp (f(x)/T)\, hi = dh/dk, 
a n (x) = (2n/ir)(n 2 + 4x 2 ) -1 , and A nm (x) = S nm S(x) + 

(1 — &nrn)[ a \n-m\ + a n+rn + 2 Ylk=l ' ^ a |ri-m| +2fe ( a; )] ■ 

The thermodynamic potentials of the host superconduc- 
tor, f2/j, and the impurity Qi, are found to be 



£lh E B cs 



dk 
2tt 



F[-u(k)] 



n i = 2e d -i{e d /2T). 



(9a) 
(9b) 



Therefore, the equation for the order parameter, 8{Q,h + 
fii)/<5A = 0, takes the form 

gL f rifc tanh(a;(fc)/2r) g Sth 
2 J 2tt w(jfe) 2A5A' 1 J 

where the first term is the standard BCS term, while the 
second term describes the impurity contribution. The 
low-temperature thermodynamics of the U — > oo Ander- 
son impurity embedded in a BCS superconductor is thus 
described completely by Eqs. (8)-(10). 

At finite temperatures, the thermodynamic BA equa- 
tions require a numerical analysis. However, at T = 
they are significantly simplified and can be solved analyt- 
ically. Indeed, one can show that the energies e(fc), 77(A) , 
and K n (X) are positive. Therefore, as in the normal al- 
loys, the ground state of the system contains only normal 



charge complexes and is described by a single equation 
for the function £(A) = £-(X)6{Q - X) + £+(A)0(A - Q), 
where Q is a single zero of the function £(A), £(Q) = 0, 

poo 

£(A) = d(A) + / dX'R(X- A')£+(A'). (lla) 
Jq 

Here, R(x) = j(dw/2-K) exp (— iux)/{l + exp \to\), and 

/oo 
dX'R(X- X')co n (X'). (lib) 
-00 

Eqs. (11) completely determine the physical properties 
of the ground state of the system, but for our purposes 
it is more convenient to derive equations describing the 
ground state directly from BAE. Introducing the "par- 
ticle" and "hole" densities of normal charge complexes, 
er(A) = 0, A > Q, and <r(A) = 0, A < Q, respectively, we 
find in the continuous limit of BAE 



1 dk(X) 1 



+ -a 2 (X - e d /2T) = / dX'a 2 {X- X')a(X') 
L J -00 

+a(X) + a(X), (12) 



2tt dX 



where fc(A) = fc( + )(A) + fc(~)(A) is the momentum of com- 
plexes. The functions a(X) and a(X) are divided into the 
host and impurity parts, i. e., cr(A) = <Th(X) + L~ 1 ai(X), 
<t(A) = (Jh{X) + L~ 1 <Ji(X). The gap-dependent part of 
the impurity energy is then given by 

^(A) = ^$, <P = J Q dX<j>(X)<j t (X), (13a) 

where 0(A) = (A — e d /2T)/(X 2 + 1/4), while the function 
<7i(X) is determined by the same equation as in the normal 
state, 

o-<(A) + cri(A) = R(X - e d /2T) + / dX'R{X - A')^(A'). 

(13b) 



At T = 0, Eq. (10) reads 
^ D dk 



1 



= gL r 

~ 2tt J 



VA 2 + fc 2 2A 5 A ' 



(14a) 



where ujd is the Debye frequency. By inserting the well- 
known solution of the Wiener-Hopf equation (13b) into 
Eq. (13a), we finally obtain 



E, 



A 2 f c 



duj 4>(—oj) exp (iuQ) 



f 

J — i 



2tt G(-)(o 
doj' R(u')G^{u') 
2m u' — uo + iO 



exp(iv'e d /2r). (14b) 



Here, <t>{oj) and R(cu) are the Fourier images of the func- 
tions 4>(X) and -R(A), and e* d = e d — 2TQ is the renormal- 
ized impurity level. The functions G^ + '(w) = G^~\— uj), 
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G(-\2mj) = V2^[(jw + 0)/e] l "/r(l/2 + iw) are analyti- 
cal functions in the upper (+) and lower (— ) half-planes. 

Eqs. (14) explicitly determine the order parameter at 
T = 0. Neglecting a small gap correction to the "normal" 
value of Q = -(1/2tt) In (pW/T), one easily finds the 
solution of Eq. (14a) in the multiimpurity case 

A = A exp (-/Zimp). (15) 

The parameter /ii mp — (cje^/r) <E>, which is proportional 
to the impurity concentration Cj, describes the impurity 
contribution. 

The qualitative behavior of the impurity contribution 
is clearly seen without a detailed study of the integral in 
Eq. (14b). Let us first consider the Kondo limit of the 
problem, where the renormalized impurity level lies much 
below the Fermi energy, — e d /2T ^ 1. The function <fi(\) 
is then positive on an essential interval of integration in 
Eq. (13a), Cd/2T < X < Q, and hence the parameter /Lti mp 
is also positive. The asymptotic estimate of the integral 
(14b) gives /ii mp ~ CiEp/\cd\- Thus, in the Kondo limit, 
magnetic impurities depress superconductivity. 

If the impurity level is shifted to the vicinity of the 
mixed-valence regime, |ej|/2r < 1, the parameter /ij mp 
changes the sign at some point e<j = id < Q ■ At > id, 
/iimp is negative, and hence the Anderson impurities, 
which in the mixed- valence regime play the role of a non- 
magnetic resonance energy level rather than that of a 
local magnetic moment, enhance superconductivity. 

In summary, making use of the BA approach, we 
have derived the exact equations describing the low- 
temperature thermodynamics of the model (2). We also 
have derived an equation for the order parameter A, min- 
imizing the thermodynamic potential of the system with 
respect to A [H. Finally, at T = we have evaluated 
the impurity part of the total energy of the system, and 
thus found an exact zero-temperature expression for the 
order parameter. 

The results obtained have a clear physical meaning. 
The ground state of a normal magnetic alloy is well 
known to be composed of the charge complexes Q . The 
appearance of a superconducting energy gap results in 
a gap correction to the energy of these complexes, and 
hence to the impurity part of the total energy of the 
system. The sign of the impurity contribution to an 
energy balance, 5Ei/SA, is different in the Kondo and 
mixed-valence regimes, leading respectively to either a 
depression or an enhancement of superconductivity. Eq. 
(15) shows that at T = the system remains in the su- 
perconducting state at any impurity concentration. Due 
to the well-developed Kondo screening of a local mag- 
netic moment of impurities, the system exhibits a Cooper 
pairs weakening rather than a pairs breaking predicted 
by the AG theory M, not accounting for the Kondo ef- 
fect. Moreover, in the mixed-valence regime, an impurity 
contribution to the density of states of the system near 



the Fermi level even dominates over a Cooper pairs weak- 
ening, and superconductivity is enhanced. 

Finally, it should be emphasized that the results ob- 
tained do not contradict a concept of gapless supercon- 
ductivity suggested by Abrikosov and Gor'kov p|. At 
T = an energy gap in the spectrum of the system may 
vanish at some critical impurity concentration, while an 
order parameter along with the parameter A remains fi- 
nite. The BA technique admits an analytical computa- 
tion of an energy gap at T = 0, that can clarify this very 
important and interesting question. 

I would like to thank A. A. Abrikosov, S. John, and 
especially P. B. Wiegmann for stimulating discussions. 



* On leave from Landau Institute for Theoretical Physics; 

current electronic address: rupasov@physics.utoronto.ca 
[1] A. A. Abrikosov and L. P. Gor'kov, Zh. Eksp. Teor. Fiz. 

39, 1781 (1960) [Sov. Phys. JETP 12, 1243 (1961)]. 
[2] For a review of earlier publications, see E. Miiller- 

Hartmann, in Magnetism Vol. 5, edited by H. Suhl (Aca- 
demic, New York, 1973). 
[3] P. Schlottmann, J. Low Temp. Phys. 47, 27 (1982); K. 

Satori, H. Shiba, O. Sakai and Y. Shimizu, J. Phys. Soc. 

Jpn. 61, 3239 (1992); 62, 3181 (1993); L. S. Borkowski and 

P. J. Hirschfeld, J. Low Temp. Phys. 96, 185 (1994); W. 

Chung and M. Jarrell, Phys. Rev. Lett. 77, 3621 (1996). 
[4] A. M. Tsvelick and P. B. Wiegmann, Adv. Phys. 32, 453 

(1983); N. Andrei, K. Furuya, and J. H. Lowenstein, Rev. 

Mod. Phys. 55, 331 (1983); P. Schlottmann, Phys. Rep. 

181, 1 (1989). 

[5] V. I. Rupasov and M. Singh, J. Phys. A 29, L205 (1996); 
Phys. Rev. Lett. 77, 338 (1996); Phys. Rev. A 54, 3614 
(1996). 

[6] V. I. Rupasov, Phys. Lett. A 237, 80 (1997). 
[7] D. Withoff and E. Fradkin, Phys. Rev. Lett. 64, 1835 
(1989). 

[8] In a pure BCS superconductor, the parameter A is both 
an energy gap and an order parameter. As in Refs. JlJ^|3j, 
we use the term "order parameter" , because A vanishes in 
the normal state of the system. In our case of a strongly 
correlated electron system, the parameter A is auxiliary. 
A real order parameter of the superconducting transition 
is only proportional to A, and may differ from both an 
energy gap in the spectrum of the system Jl]] and from the 
parameter A. 

[9] See, for example, M. Tinkham, Introduction to Supercon- 
ductivity (McGraw-Hill, New York, 1996). 

[10] P. B. Wiegmann, Phys. Lett. A 80, 163 (1980). 

[11] H. Shiba, Prog. Theor. Phys. 40, 435 (1968). 

[12] Eqs. (10) and (14a) are not related to integrability of the 
model. They are obviously valid for any dilute magnetic al- 
loy, provided superconductivity of a host metal is described 
by the BCS theory. 



4 



